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Precise measurements of polarized electro-production and 
Drell-Yan in the deep inelastic limit will soon provide first 
information on the higher twist parton distributions gr(x) 
and hi,{x). Sum rules for higher- twist structure functions 
are only valid provided the corresponding Compton ampli- 
tudes satisfy un-subtracted dispersion relations. Subtracted 
dispersion relations have to be used when the (real part of 
the) forward scattering amplitudes does not fall off rapidly 
enough for v — > oo (fixed Q 2 ). Formally, such subtractions 
lead to ^-functions at the origin in the parton distributions, 
which are not accessible to experiment, and the integral over 
the data fails to satisfy the sum rule. The (5-functions in the 
parton distributions can be identified with the zero-modes 
that appear in light-front quantization. An explicit infinite 
momentum boost identifies these soft quark modes with low 
momentum contributions to the self-energy. 



I. INTRODUCTION 

With the advent of precise deep inelastic scattering 
(DIS) experiments it will soon become possible to mea- 
sure the proton's polarized higher-twist structure func- 
tions gr{x) and }il(x). 9t(x) = g2{x) — g\{x) is obtained 
in polarized electro-production [jp, where one finds for 
the parallel-anti-parallel asymmetry 



da" 



Ana 



dq 2 dE' dq 2 dE' E 2 q 2 



[{E + E' cos 9) MG X + q 2 G 2 ] , 
(1.1) 



where E, E' are the lab energy of the initial/final lep- 
ton, f, J, denote nucleon/lepton felicities and 9 is the 
lab angle of the final lepton. In the parton model, 
vG\ scales while only v 2 G2 has a finite scaling limit as 



Q 



oo, v = E — E' — > oo and xbj 



l /2Mv 



fixed (the Bjorken limit). Using the operator product ex- 
pansion one finds (modulo logarithmic scaling violations) 



W'vGx ^ E,eJfl?(x BJ ) + 9l(-x Bj ) and Mv 2 G 2 
J2q e "q92( x Bj) + sfC - X -Bj){] where the parton distribu- 
tions g\ and g\ are defined through the light cone corre- 
lations B 



Bj 



(1.2) 



2 [g q (x)p»(S ■ n) + g q T (x)S£ + M 2 g q (x)n»(S ■ n)] , 



where n 2 = p 2 = 0, n + = p = 0, P^ = p^ + ^-n**, 
=p»(S-n) + nV{S-p) + S£ and g q T (x) = g\{x)+gl{x). 
The different scaling behavior of G\ and G2 as p + —* 
(Breit frame) is reflected in the coefficients in Eq.( |l.2[ ), 
where the coefficient of g\ is inde pendent of p + , while the 
contribution from g\ in Eq.(1.2) grows linearly with p + . 
This is typical for operators of different twist. Clearly, 
the non-leading role played by G2 makes it rather diffi- 
cult to separate them from the leading twist term G\ (and 
from 1/Q 2 corrections to G\). The HERMES experiment 
at DESY will be the first attempt for a precise determi- 
nation of g^{x) ||, which has received special attention 
due to the existence of the Burkhardt-Cottingham sum 
rule 1 



dxg q (x) = / dx [g q (x) + gl(x)] = 0, (1.3) 



which seems to follow easily from the OPE result ( |l.2| ) 
HP!- To see this, one can integrate Eq.( |1.2j ) over x [using 
the fact that the r.h.s. has support only for a; € (—1; 1)], 
where we find 

(PS| g -(0)y 7 59(0)|PS) = (1.4) 
2/lr [9 q i(x)p"(S ■ n) + g q T (x)S£ + M 2 g q (x)n»(S ■ n)] . 

fro m ro tational invariance it follows that the l.h.s. in 
Eq.(L4) is proportional to the spin vector 5^ and thus 



the g q (x) must satisfy 



dxg^(x) 



1 r i 

dxg\{x) - 
-1 J-i 
1 ,1 

dxg\(x) = 2 dxgl{x). 
-1 J-i 



(1.5) 



The first of these conditions (1.5) is abovementioncd 
Burkhardt-Cottingham sum-rule, while the second is a 
(probably useless) sum-rule for 33 .^] 



2 Useless, because no experiment is known that could access 
l gl{-x) is sometimes denoted gf (as) in the literature. ff3 with sufficient precision and because Regge arguments sug- 

gest that the sum rule is highly divergent. 
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The way I presented the derivation of these sum rules, 
it seems that they are mere consequences of rotational 
invariance and hence there should be no question about 
their validity in QCD. However, there are several sub- 
tleties if one wants to turn these sum rules for the g^s 
(the light-cone correlators) into sum rules for the GVs 
(the experimentally measured cross sections). Discussing 
these subtleties will be the main purpose of the rest of 
this paper. Here I only mention that the OPE result is 
derived using dispersion relations for the lepton-nucleon 
forward scattering amplitude and in general there is al- 
ways the possibility that the dispersion relation does not 
converge or that the dispersion relation requires a finite 
subtraction. 

Similar results apply for the chirally odd parton distri- 
butions hf(x) defined as M 



(1.6) 



^e iXx (PS\q(Q)a^i l5 q(Xn)\PS) 
Ztt 

= 2 [h\(x) {S^p v - S$.ff) /M 
+h g L (x)M {p»n v - p v n^) (S ■ n) 
+ h q 3 (x)M(S£n u - S£n")], 

where h q L {x) — h\(x) + h\(x)/2. Since they are chi- 
rally odd, they cannot be measured in inclusive electro- 
production. One promising way of measuring them is the 
polarized Drell-Yan process Q (see also Ref. Q): If both 
target and beam are longitudinally polarized, one finds 
for the spin assymetry || 



A LL = 



E q e 2 q gl(x)gl(-y) 



(1.7) 



If both target and beam are transversally polarized, || 
_ sin 2 flcos20E g effi(a:)/if(-i/) , , 

and if one is longitudinal and the other transverse, PJlO| 
2 sin 20 cos M 



A 



LT 



I + COS 2 8 ^/Q2 

E q 4 (gf(x)ygU-y)-xhl(x)hU-y)) 



(1.9) 



E q e 2 q f?(x)f!(-y) 



(1.10) 



An alternative way of measuring h\ is through semi- 
inclusive polarized electro-production 

Similarly to the sum rules for the o| one can derive 



sum rules for the h\ by integrating Eq.(1.6) and one finds 
l&BJ, using again rotational invariance, 



i ,i 

dxh\{x) — I dxh q L (x) 
l J-i 
*l pi 

dxh\{x) = 2 1 dxh\{x). 



(1.11) 



As is the case for the 33-sum rule, the /13-sum rule is most 
certainly useless, but since measurements for h\ and hi, 
are u nder way, it will soon be possible to test the /i^-sum 
rule ( 1.11 ) experimentally. In principle, the same com- 
ments that I made above for the gx sum rule also apply 
for the hh sum rule. However, as I will discuss in Sec- 
tion III, it seems to be violated already in perturbation 



theory. 

The paper is organized as follows. In Section ||, I 
will explain the connection between the light-cone corre- 
lation functions and the experimentally measured quan- 
tities and why the formal sum rules — derived for the 
light-cone correlation functions — may be violated when 
app lied to the measured structure functions. In Section 
III , I will show that the h% sum rule is violated in QCD. 
In the following sections I will study some model filed 
theories and simple (spin-independent) higher-twist par- 
ton distributions and I will attempt to shed some light 
on the physics of such violations of sum rules. In the 
Appendix, I will draw connections to the notorious zero- 
modes plaguing the light-front community. 



II. OPERATOR PRODUCT EXPANSION, 
DISPERSION RELATIONS AND ALL THAT 

Deep inelastic scattering (DIS) experiments are per- 
formed in the region < Q 2 < 2Mv [O]. The operator 
product expansion is an expansion around Q 2 — 00 and 
must diverge for Q 2 < 2Mv because of the singularities 
of the hadronic tensor (the physical discontinuity along 
the cut). The consequences of these statements are eas- 
ily demonstrated by considering some hadronic tensor 
T(x,Q 2 ), where x = Q 2 /2Mv. To keep the discussion 
as general as possible, we will consider here some generic 
forward scattering amplitude, without specifying the cur- 
rent. This will allow us later to draw very general conclu- 
sions. A very similar discussion, for the specific example 
of gi{x), can be found in Ref. [Q. 

As an function of x, T(x, Q^is analytic in the com- 
plex plane, cut along < x < 1. The discontinuity of 
the imaginary part along the cut is (up to trivial kine- 
matic factors) the experimentally measured cross section, 
which we will denote by G(x, Q 2 ). In order to apply the 
OPE, one needs to relate T(x, Q 2 ) for |x| > 1 to T(x, Q 2 ) 
for I a; I < 1. From the theory of complex functions it is 
well known that an analytic function is determined by its 
singularities in the complex plane — up to a polynomial! 
This implies (|x| > 1) 



T(x,Q 2 



x 



dx'- 



^G(x'), (2.1) 



where p Q 2 ) is a polynomial in - whose coefficients do 
in general depend on Q 2 . In the language of Regge the- 
ory, such polynomial subtractions are called "fixed poles" 
p3|. Note that we specialized here on an odd function 



2 



of 1/x, the generalization to the even case follows anal- 
ogously. Usually, only a few coefficients of the polyno- 
mial are allowed to be nonzero, since otherwise unitarity 
bounds are violated, but typically this places no restric- 
tion on the lowest few terms. 

For x > 1 one can apply the OPE and one finds 



T(x,Q 2 



E 



i 



(2.2) 



n=l,3,5,. 



where the coefficients are related to matrix elements of 
local operators. If one allows for generalized functions, it 
is always possible to find a function g(x'), such that 



dx'x' n 1 g{x) Vn = l,3, . 



(2.3) 



In the case of the moments appearing in the OPE for 
DIS this function typically is a Fourier transform of a 
correlation function along the lig ht-cone, such as the gi 
and hi introduced in Eqs.(1.2,1.6). 

In order to demonstrate what happened if the sub- 
traction polynomial does not vanish, let us study an ex- 
ample, where the polynomial is finite^ and fir st o rder^ 
P (x>Q 2 ) = ^Pi(Q 2 )- Upon inserting p in Eq.(|l]) and 
introducing the moments of the measured structure func- 
tion 



b n = 



1 



, / /n-l 

ax x 



G(x') 



(2.4) 



j(x) + g(-x) = Pl S{x) + -G(x) (2.6) 

7T 



A few more comments are in order here. So far, I have 
only discussed what happens if there is such a subtrac- 
tion. In the rest of the paper I will show (studying QCD 
perturbation theory as well as some toy models) that this 
is actually a very common phenomenon. In these exam- 
ples I will furthermore demonstrate, that the light-cone 
correlations indeed contain 5-functions, even when the 
"data" is very smooth, and that the coefficient of the 
(5-functions indeed coincides with subtraction constants 
in dispersion relations. Furthermore, I should empha- 
size that the ^-functions do not arise from some obscure 
Q 2 — > oo limit but from inverting a moment transforma- 
tion at fixed and finite Q 2 . 



III. PERTURB ATIVE ANALYSIS OF H L AND G T 

So far, the discussion has been deliberately abstract. 
Even though we have seen what happens if the dispersion 
relations contain subtractions, we have not yet shown 
whether this actually happens. As a first step in this 
direction we will perform a perturbative analysis of h^ (x) 
and gr(x) in this section. Let us first look at the 0(g 2 ) 
corrections to hh{x) for a target that consists of a quark 
of momentum P and spin S, dressed by one gluon loop 
(Fig#. 



one thus finds 



an — b n tl — 3,5,. 
ai = pi + h . 



(2.5) 



Such a result, if p\ is nonzero, has a two major conse- 
quences: 

• Any sum rule, derived for a\, fails to be satisfied 
for b\. Whether or not the sum rule derived us- 
ing the OPE applies also to the experimentally ob- 
served structure functions depends on whether or 
not the amplitudes satisfy un-subtracted dispersion 
relation. 

• Even when G(x') is a smooth function, since all 
moments except the lowest moment of G(x) /ir and 
g{x) + g(—x) are the same, it must be that g{x) 
contains a (^-function at the origin (which obviously 
afflicts only the lowest moment) 



3 N ote that in general it may happen that the integral in 
Eq^^l) does not converge and hence one needs an infinite 



subtraction. Such a case may be considered as a limit of the 
situation assumed here. 

4 Typically, higher orders are not allowed anyway because of 
unitarity constraints. 



FIG. 1. One loop correction to gr{x) and Iil(x) in QCD. 
The cross represents the insertion of the light-cone correlator. 

For simplicity, we will assume P±_ = 0, but the gen- 
eralization to nonzero trans vers e momenta is straightfor- 
ward. From the definition (1.6) one finds 
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hL{x)S ~M 



-g 2 u(P,S) 



d 4 k 



S(k + - xp + )Y' 



where 



«75" 



• IE 



'fu(P,S)D ia ,(P-k), 



(3.1) 



DM 



■ IE 



nq 



(3.2) 



is the gluon propagator i n th e light-front gauge 5 and m is the current quark mass. The term S(k + — xp + ) arises from 
the A-integration in Eq.(.L6) which projects out quarks of momentum fraction x. 

Straightforward application of the 7-matrix algebra yields (using Pu(P, S) = m u{P, S) and dropping terms odd under 
k± -> -k±) 



h L (x) 



~M 



g 2 u(P, S) 



d 4 k S(k + — xp + ) (—4m(k + -f —k r y + )ij5 



(2tt) 4 - k) 2 + is 



1 [27717+175 + 2k + <r + 175 



(fc 2 



Most of the terms in Eq.(3.3) are harmless as k + — > 0. 
The only troublesome piece arises from the term propor- 
tional to k~ . In order to study its contribution further, 
one can make use of the algebraic identity 



P- 



2{P^ 



(P-k) 2 

2{p+-k+y 



(3.4) 



One can easily verify that the first two terms on r.h.s. of 
Eq.(3.4) are well behaved as k + — > 0, when inserted in 
Eq.(3.3). The only troublemaker is the term proportional 
to (P — k) 2 in the numerator. Using 



dk- 



1 



i 6(k^ 



2ir {k 2 — m 2 + ie) 2 
one finds for the singular piece 

h s i na {x)MS + = g 2 mS+P^ 



2k 2 ± 



(3.5) 



2 mS+ 
9 o{x) 



dk+d 2 k ± S(k+ - xP+)S(k+) 
(2tt) 3 kl+m 2 
! 2 fc ± 1 



(2tt) 2 k 2 



2 mS+ A 2 ± 
g -r —6{x)\o<g — 
vk 1 m z 



where is some cutoff and where we used S ■ P = 0, 
i.e. S~ = -S+M 2 /2P+. The factor S+ on the r.h.s. 
arises from the matrix element of «7 + 75- The precise 
numerical result is not so important here. What is im- 
portant is that straightforward evaluation of our pertur- 
bative model yields a term proportional to a 6 function 
for tiL(x). 



IE) 



k+ 



IE 



u(P,S). 
(3.3) 



The singular term derived here is explicitly proportional 
to the current quark mass, i.e. only for strange or heav- 
ier quarks one expects a sizable contribution, and thus 
violation of the /12-sum rule. However, it is not known 
whether this is only a lowest order artifact or whether 
this is a more general result. Certainly, using hard-soft 
factorization, this result implies a S(x) term with nonzero 
coefficient also in the full (i.e. nonperturbative) result for 
Jil{x) and thus the /i^-sum rule should be violated for 
strange and heavier quarks. 

It should be emphasized that while we have identified 
the quark mass as one source of violation for the /i^-sum 
rule, it is by no means clear whether or not there are are 
other sources of further violation which would also affect 
light u and d quarks. 

Before we are going to investigate the physical origin 
of the (5-function, let us first take a look at QCD pertur- 
bative corrections to gr(x) for a quark. 



(3.6) 



4 



gr(x)Si 



-g 2 u{P,S) 



d A k 



(2tt) 



IE 



"7 «75" 



From the purely formal point-of-view the (^-function in 
Iil (x) arises because of the term proportional to the light- 
front energy fc~ = (k° — k 3 )/\/2 in the numerator, which 
cancels the energy dominator for the gluon propagator. 
Armed with this insight, we can immediately focus on 
the potentially dangerous terms in gx arising from the 
k~ in the numerator of the fermion propag ator| 

(3.8) 



7ai r y + j x ij5fn + mj x i-f 5 k 7 + ] 7^ = 0, 



i.e. because of the Dirac algebra, the coefficient of the 
potentially singular term in gr{x) vanishes. In contrast 
to the result for hi,(x), at least in one loop order, radia- 
tive QCD corrections do not give rise to (5-functions in 
9t{x). 



-^u{P,S)D^{P-k). (3.7) 



IE 



we will consider only the symmetric part of the "H — " 
com pone nt. A typical Feynman diagram contributing to 
Eq. (U) to <D{N°) (leading order) is depicted in Fig. §. 



IV. LONG-RANGE CORRELATIONS ALONG 
THE LIGHT-CONE IN A SIMPLE TOY MODEL 

In the previous section, we have seen that sum rules 
for higher-twist parton distributions can be easily vio- 
lated. However, so far, we lack an intuitive understand- 
ing of the effect. In this (and the following) section we 
will turn our attention to simple toy models, where the 
same effect happens, but where it is easier to analyse. 
For the same reason, we will study the (spin indepen- 
dent) scalar higher- twist distribution e(x) in the rest of 
the paper. Both, working in 1+1 dimensions and study- 
ing scalar distributions will simplify the algebra involved 
considerably and thus allow us to shed some light on the 
essential physics. The general conclusions should not be 
affected by these simplifications. One of the most sim- 
ple examples for violations of sum rules for higher-twist 
distributions is a(l + l)-dimensional model of the Gross- 
Neveu type HJTq] 



C = Ui(i $ — too) Ui + Si (i 



mo) Si - —UiUiSjSj 

N r 



(4.1) 



where i,j — l,...,Nc are "color" indices. Here we will 
restrict ourselves to leading order in l/N c . Suppose we 
probe a it-quark with an external vector current which 
couples only to strange quarks. From the hadronic tensor 



d 2 xe zqx (u,p\T (s^s(x)sYs(0))\ u,p) 

(4.2) 



FIG. 2. Leading order l/Nc contribution to T£J U in the 
Gross-Neveu model. The dotted lines indicate the flow of 
momentum through the external currents S7 M S. 



The OPE yields for Q 2 



-q- — > 00 



T + r + t -+ 

t S/U S/U 



2M 



2 1 dx'- 

1 : 



2 c s/u\ 



(4.3) 



where 

/oo 
dXe^ x (u,p\s{0)s(Xn)\u,p) . (4.4) 
-00 

Here M is the constituent quark mass JlJ,|l^]. Clearly, 
e.i/n(x) satisfies a "sum rule" (the sigma-term sum rule 

£3) 



dxe s/u (x) = —(u\ss\u) 



(4.5) 



which is non-zero and0(Ar?) @. However, any attempt 
to verify Eq. (4.1:) by studying physical cross sections is 
doomed to fail, since STf^ = to D(N°).f\ 

These seemingly contradictory results are consistent 
with di sper sion theory. The crucial point is that the OPE 
result ( |4.3| ) is derived in a kinematic regime where T^ v 
is purely real |a;| > 1 and which is not accessible in DIS. 
A relation between T for \x\ > 1 (where T is purely real) 
and for |a;| < 1 is established by means of a dispersion 
relation (x > 1) 



One can easily verify that terms with k~ in the numerator 
from the gauge field propagator are accompanied here by an 
explicit factor tx k + in the fermion propagator and are thus 
harmless. 



7 This can be easily verified by applying the cutting rules to 
Fig. 1. Since q 2 < 0, at least one of the cut lines will have a 
negative invariant mass. 
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T 8/u (x)=p(-)+ - [ dx' X l2 %T s/u {x') , (4.6) 

\X J 7T Jq X X 

where p(l/x) is an odd polynomial in 1/x which cannot 
be determined from dispersion theory — it corresponds 
to subtraction constants. If p would vanish, then a com- 
parison between Eqs. (4.6) and (|4.3|) shows 



~ST s/u (x') = 2M 2 [e(x')+e(-x')] . (4.7) 

7T 

However, if for example, p(l/x) — c/x then 
-5sT s/u (x') = c-S{x') + 2M 2 [e(x') + e{-x')} . (4.8) 

Since 3T, which is measured experimentally, does not 
contain ^-functions, we expect them to be contained in 
e(x') +e(— x'). The latter will now be determined explic- 
itly: obviously (Fig. ||) 



(u,p\s(x)s(y)\u,p) 
d 2 k 



(4.9) 



9 



(27T) 



-u{p)u(p)tr 



-iky 



m + ie ft — to + ie 



i.e. 



s/uW)=9 2 / 77" 



dX 



:k + 



,i\x 



2tt 



dk+dk- (k 2 +m 2 )e *F 



(2tt) 2 (2k+k- - to 2 + lef 



(4.10) 



FI G. 3. Leading order 1/N C Feynman diagram for 
Eq. Q. 



Since the exponent in Eq. (4.10) does not contain fc~, 
one can always close the contour in the complex fc _ -plane 
such that no pole from the energy denominators is en- 
closed. For k + 7^ the surface term (from the semi-circle 
in the complex k~ -plane) vanishes and one finds 



dk~ 



(2k+k- - m 2 + ieY 



oc 



5{k+) 



(4.11) 



[The coefficient of proportionality is non-zero and cal- 
culable. It turns out to be such that e s / u {x) satisfies 
Eq. (4.5), as one can verify by comparison]. Therefore 



u (x') oc 5{x') 



(4.12) 



Mathematically the 8- function in Eq.( [4.12 ) arises as a 
consequence of long-range correlations along the light- 
cone. In the Gross-Neveu model, to leading order in 
1/N C , the light-cone correlator (u,p \s(0)s(x~)\ u,p) is in- 
dependent of x~ — which is a consequence of Eq. (2.11). 
Physically, this is related to the scalar coupling in the 
i-channel. This remark becomes more clear if we define 
the light-cone distributions by means of an infinite mo- 
mentum boost fllTf . 

Consider the equal time correlator (p = p 1 , k = k 1 ) 



P P (k') 



dx\ 
~2w' 



ikxi 



u,p\s(0)s(xi)\u,p) 



(4.13) 



which measures nothing but the conventional (i.e. not 
light-cone) momentum distribution flowing through a 
strange current insertion (with scalar quantum numbers). 
It is related to the light-cone distribution through 



e(x) = lim p ■ p p (x • p) 

p — >oo 



(4.14) 



Now, as a consequence of the scalar t-channel coupling, 
p p (k') is actually independent of p. One finds 



p p (k) 



dk Cl 
~2^ 



-tr 



k 2 - 



k 



(k 2 + to 2 ) 



3/2 



(fc 2 + M 2 ) 



?,/2 



'to -> M" 

(4.15) 



where we subtracted the contribution from a heavy reg- 
ulator, in order to rende r (u\ss\u) = j dk p p (k) finite. 
Clearly, p p (k) (Eq. fl4.15| )) is a localized function with fi- 
nite area. Measuring all momenta in units of p (wh ich i s 
after all what the area preserving mapping in Eq. (4.14) 
does) and performing the limit p — > oo thus gives rise to 
a distribution of zero width but finite area — a represen- 
tation for a (5-functionJ^| 

The fact that p p (k) is independent of p, to leading or- 
der in 1/N is quite obvious, since the leading order in the 
1 /N approximation is similar to a mean field approxima- 
tion. The expectation value of a scalar density, taken in 
a plane wave state, is a constant — independent of the 
momentum. Therefore, the u-quark gives rise to a source 
for s -qu arks [through the point-like scalar coupling in 
Eq. (4.2 )1 which also does not depend on its momentum 
(ET5). 



Thus one might suspect that the unique properties of 
the interaction term in Eq, (^^) are the reason for our 
"weird" results: 



Using various test functions can can verify that Eq. (4.14) 
yields a "pure" 5-function and is not contaminated with 
higher multipoles, like 8" . 
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• smooth (actually identically zero in this exam- 
ple); 

• 3T violates the sum rule (the sum rule converges 
without problems); 

• e(x) contains a S- function; and 

• e(x) satisfies the sum rule. 

In particular, one might suspect that "more realistic" 
field theories (like QCD3 + i) behave differently. In the 
rest of this article we will demonstrate that some of these 
properties are actually quite generic for higher-twist dis- 
tributions in many field theories (also QCD 3+ i). 

A non-trivial and non-perturbative example has al- 
ready been given in Ref. [13], where the twist-4 parton 
distribution was analyzed in the context of QCD1+1 
(N c — > oo). There, QT~~ (which scales against for 
non-zero x in the Bjorken limit) violates the sum rule 
which is derived for f±. Furthermore, fi contains a 5- 
function, while lim QT has only a mild (integrable) 

singularity for x — > 0. 

In the above discussion, the presence of ultraviolet di- 
vergences in the Gross-Neveu model forced us to intro- 
duce a regulator (cutoff). In order to keep the discussion 
simple the issue of renormalization (and elimination of 
the cutoff dependence) has been avoided so far. This will 
be done in the following. To leading order in 1/N C one 



still finds 3T, 



In 



amplitude one finds H], (Fig.@) 



0. for the real part of the renormalized 



9 , x f d 2 k 



Mp 



1 



#+ i-M F 



(4.16) 
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where Mp is the physical fermion mass. g aq q is the phys- 
ical quark-cr- meson coupling and D& (0) is the a- meson 
propagator evaluated at p 2 = (note that no momen- 
tum flows through the tr-meson, since we are considering 
a forward amplitude). Therefore, 



T, 



=v(t+- u+ t;+) 



M> f 1 , (Q 2 ) x(l-x 
c—-7r- I ax 



o [Q 2 x(l -x) + M 2 f 



vM 2 (<?_ 
Q 2 \Ml I ' 



(4.17) 



where c is some numerical constant. If one now intro- 
duces a Q 2 -dependent parton distribution eg/u(x' 7 Q 2 ) 
via the moment expansion 

1 f 1 

T s/u (x,Q 2 ) = 2M 2 — dx'x' 2 e s/u (x',Q 2 ) 

,,-n OA X J -I 



one finds the unique solution for the even part of 
e s /u{x',Q 2 ) 



v 



u y& i 



■'s/u 



4M f 



2M F 



(Q 

M 

(4.19) 



Since T^ v is independent of v, for all Q 2 , one has 
v ■ T cx l/x(x — Q 2 /2Mfv), Therefore, the (^-function is 
present for all values of Q 2 . The coefficient multiplying 
S(x') depends on log Q 2 /M F which reflects the running 
of the coupling constant in the Gross-Neveu model. 



V. BOOSTING SLOWLY TO INFINITE 
MOMENTUM 

In this Section we will continue our analysis of higher- 
twist distributions in various field theories, although we 
will now use perturbation theory to the lowest non-trivial 
order only. Of course, parton distributions in QCD can- 
not be calculated perturbatively, but this is not what we 
attempt to do here. The main motivation for this pertur- 
bative study is to shed some light on possible mechanisms 
for the violation of un-subtracted sum rules. The use of 
perturbation theory will allow us to explicitly compute 
various scattering amplitudes and correlations analyti- 
cally. This is very important if one wants to understand 
the reason for the failure of the sum rules mathematically 
as well as physically. 

The example we will be studying at is again a 1 + 1- 
dimensional model: massive fermions coupled to massive 
bosons with a chirally invariant Yukawa coupling (linear 
tr-modcl) 

C = ijj (i — m) tp — a (□ + m 2 ) a — n (□ + m 2 ) 7r 

+ ■ytp (a + ij 5 n) tfj , (5.1) 

where, for simplicity, all masses are taken to be equal. 
The flavor structure [which is suppressed in Eq. (5.1)] is 
assumed to be such that only Fig. || contributes to order 



,2 to the hadronic tensor 



T^ = i / d 2 xe lqx (u,p\T {d^siO^dix)}] u,p) 



(5.2) 

This can, for example, be accomplished by having only 
the flavor combination ud and du couple to the mesons. 
Straightforward application of Feynman rules yields 



(4.18) 
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one finds 



—Aim 2 ^ 2 u{p) 
16m 2 7 2 



d 2 k 



(2tt) 2 ( fc 2 _ m 2 + ie f f^ k + q f _ m 2 + ( p _ k f _ m 2 + 



■u[p) 



4tt 



dy 



dz- 



(zm 2 — yp ■ q) (1 — y — z) 



[y(l — y)q 2 + z(l — z)p 2 + 2yzp ■ q — m 2 + ie 



For the application of dispersion relations it is helpful to perform an integration by parts in Eq. (|5.3|) 



2 ? 

m 7 



dy Q-JZ 

o [y(l - y)q 2 - m 2 



9 ? 

m 7 



dy 



1-2/ 



dz 



/ 2m 2 (l - y - 

\ 1)3 



1 

T3 2 



(5.3) 



(5.4) 



where 

D = y(l — y)g 2 + z(l — z)p 2 + 2yzp ■ q — m 2 + ie 



(5-5) 



2me d/u (x) 



m 2 7 2 ' x 



1 - ) (1 - x)d(x)d(l - z) 



47T 

2m7 2 <5(ir) 

47T TO 2 



[ra 2 cc(l 



7 2 J x 2 6>(a;)6>(l - x) 
2^i 1 [x(l - x) - if] 



5{x) 



(5.8) 



FIG. 4. Box diagr am contribution to T M " (5.2) in the linear 
cr -model [Eq. (£jj)]. 



where the <5- function came from the last term in Eq. (fTq). 
It should be emphasized here that the 5- function in e(x) 
is closely related to the non-covariant counterterms in 
the Hamiltonian formulation of light-cone field theories 
]l8|-p0| . It should also be emphasized that the same trick 
works in 3 + 1-dimcnsions, i.e. there is nothing special 
about two dimensions here. 

Mathematically, the ^-function appeared here because 

IB 



Most importantly, the first term in Eq. (5.4), which is 
purely real for q 2 < 0, does not depend on p ■ q. The 
presence of this term thus ruins the application of an un- 
subtracted dispersion relation to (T -1 +T h ) /2. As 
explained in Section |n|, the presence of subtraction con- 
stants in the dispersion relation is correlated with the 
presence o f a ( 5-function in the light-cone momentum dis- 
tribution (2.6): 



dkr 



1 



2vr [k 2 - m 2 + ie)' 



= S(k+) 



1 



2m 2 



5(k 



fk 

2?r (k 2 -m 2 + ie) 2 

-) . (5.9) 



With 



2e d /„(x) 



d 2 k 



HP) JiS(^r-x)u(p) 



From the physics point of view it is very instructive to 
consider a careful boost to infinite momentum. In anal- 
ogy to Eq. ( 4.13j ) we consider the momentum distribu- 
tions of (i-quarks, in the wave function of a dressed u- 
(5 6S uar ^ °^ momen tu m p. One finds for the interaction in 
Eq. * 1 1 



(2tt) 2 (fc2 - m 2 + ie f (( p _ k y _ m 2 + ie y) p (k) 



TO7 

~2^~ 



and 



fi = k + j + k 7 + = k + j~ 



(5.7) 



(k 2 
l 

dy- 



2\3/2 

vnr) 

(k - py) 2 + m 2 (y 2 + \y 



(5.10) 



o \(k — py) 2 — y(l — y)m 2 + rri 



15/2 



2(p+-k+) 



, (p — k) — m 2 , 

7 - -77T-r — rz^7 



2(p+-k+) 



Although there is no elementary scalar exchanged in the 
i-channel, p p (k) contains a term that does not depend 
on p at all. Upon substituting k — ► xp and performing 
the infinite momentum boost, this term piles up to a 5- 
function (again one verifies that no 5" components are 
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present) as shown in Fig. ^. Actually lim pp p (xp) gives 



p — >oc 



a result identical to the light-cone calculation <\5 



FIG. 5. The momentum distribution of d-quarks in a mov- 
ing u-quark (gjj) for various values of the u-quark momentum. 
p ■ p p (x) is plotted in units of -y 2 firm. 

Although the linear cr-model does not yield a point- 
like scalar coupling in the t-channel, we observe a similar 
effect as in the GN-model: a part of the sea quark wave 
function does not follow the boosted valence component. 
In coordinate space language this means that only part 
of the sea quark distribution becomes Lorentz contracted 
when the u-quark is boosted, while another component 
looks the same in all frames! In 3+1 dimensional exam- 
ples (similar to the one above) this means that dressed 
particles do not necessarily become "pancakes" as p — ► oo 
and some component of the wavefunction retains a fi- 
nite longitudinal extension. Whether that component is 
"visible" in some experiment depends of course on the 
quantum numbers of the probe and, in a manner which 
is not yet understood completely, on the details of the 
underlying theory. 

If one calculates only 



(5.11) 



(Fig. g), one finds no indication whatsoever about the 
presence of a 5-function in e(x). For all values of Q 2 , 3T 
vanishes smoothly as XBj = Q 2 /^P • q —> 0.0 Also the Q 2 
evolution does not rise to any "peak at small x." 



9 The divergence for large x in Fig. |^ is a purely 1 + 1- 
dimensional artifact and arises from phase space factors. 



FIG. 6. p ■ (T+- + T _+ ) 1 2 (in units) as a function 
of x — Q 2 /2pq for various values of Q 2 . Note the perfectly 
smooth behavior near x = . 

Furthermore, as Q 2 — > oo, 3T approaches e(x) (Fig. 4) 
as long as x ^ 0. However, even though is very 
smooth as x — ► 0, the naive version of 



2m / dx e c i/ u (x) — (u \dd\ u) , 



(5.12) 



where one replaces 2me,j/ u (x) by fails. This should be 
obvious since (u|dd|u) is negative to order j 2 in the 
above model, while 3(T) is always non-negative (from 
unitarity). 

A different perspective is reached by computing the 
real part of T as well. In the Bjorken limit one finds 

V ' 2 



i 







[x(l -x)- I] 

'dx'^l- 
x — X 
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27T2; 



(x > 1) (5.13) 



which is consistent with the result from the OPEj^l 

,e d /u(x') 



3?T(x) = m 



dx 

! X — X' 



(5.14) 



However, the presence of the subtraction term in ( 5.13| ) 
means that one can identify ^sT(x') with e(x') only for 
x' t^O. 

It is instructive to re-analyze the situation from the 
point of view of the infinite momentum frame.^ As 



10 Notice that, since the currents in Eq. (3.2) are not charge 
conjugation eigenstates, e(x) and T(x) do not have a definite 
transformation property under x — > —x either. 

11 Of course, by making p in Eq. (3.8) large, one boosts the 
target and not the observer. However, this is equivalent. 







p — > oo one almost recovers the naive picture, where all 
virtual constituents carry a positive, finite fraction of the 
hadrons' momentum: pp(xp) drops sharply at x = 1 for 
p = 10 3 in Fig. 4 and vanishes for x < 0. This is con- 
sistent with the absence of vacuum fluctuations in the 
oo-momentum frame because this would require at least 
one particle carrying a negative momentum fraction |L7| . 
However, some component of the wave function of the 
u-quark "lacks behind" in the boost ( |5.10 ) . This can be 
summarized in the following physical picture. As p — > oo 
most of the soft modes decouple from the system, which 
is partly responsible for the simplified dynamics in the 
infinite momentum frame. Of course these modes don't 
completely disappear but get concentrated in the region 
near x = 0. While leading twist distributions couple only 
weakly to this component of the wave function, there is a 
stronger coupling of the higher-twist distributions. That 
is why leading twist distributions do not contain these 
^-functions. A very similar and related effect was first 
observed in the context of QED3 + i 21 1. There it was 



shown that certain connected vacuum graphs do survive 
the p — > oo limit. Since e{x) is, up to the momentum 
projection, proportional to a mass insertion, it should be 
clear that the non-covariant piece in the self-energy of an 
electron in QED — which is all that survives from the 
z-graph as p — > oo — is (up to mass derivative) propor- 
tional to the 5-function contribution one finds for e(x) in 
QED. 

I should emphasize again that, if such a situation ap- 
plies to experimentally measured higher twist distribu- 
tions, one would not be able to "see" any hint about the 
singular behavior of the light-cone distribution by just 
looking at the structure function. 



VI. SUMMARY 

We have investigated higher-twist distributions in a 
variety of field theories. This included perturbative ex- 
amples in 1 + 1 as well and 3 + 1 dimensions]^] In 
most cases the naive sum rules for the lowest moments of 
higher-twist distributions are violated when summed over 
the "experimentally determined"^ distributions. For 
non-zero xb), these structure functions scale towards 
the Fourier transformed quark-quark correlation in the 
hadron along the light-cone. The phenomenon which 
causes the failure of the sum rules for higher-twist struc- 



12 For a non-perturbative example in 1 + 1 dimensions, see 
Ref. jiij for a discussion of fi in the context of QCDi+i 
(N c — » oo). See also Ref. for a discussion of the scalar 
density in the sine-Gordon model. 

13 Actually, the imaginary part of some scattering ampli- 
tudes. Of course, there can only be Gedanken experiments 
in model field theories. 



ture functions is a <5-function at the origin in the light- 
cone distribution. Structure function and light-cone dis- 
tribution coincide for non-zero x, while the 5-function is 
absent in the structure function. The sum-rules are valid 
when one includes the 5- function, i.e. when applied to 
the light-cone distribution. The structure functions then 
cease to satisfy the sum rule because they do not contain 
the 5-function. In that sense the <5-function destroys the 
sum rule. From a theorist's point of view, the 5-function 
restores it, i.e. when added to the structure function the 
combined result should reproduce the sum rule. The last 
point may however be of purely academic interest, since 
one cannot measure the light-cone correlation directly. 

The relation between light-cone correlations and struc- 
ture functions is usually based on the operator product 
expansion combined with dispersion relations. Unless 
one has a theoretical or experimental constraints on the 
high-energy behavior of the real part of the scattering 
amplitude, subtraction constants cannot be ruled out. 
In fact, the coefficient of the 5-function is proportional 
to such a high-energy subtraction. This also shows that 
the long-range correlation along the light-cone are con- 
sistent with the OPE. 

In coordinate space language this means that only part 
of the sea quark distribution gets Lorentz contracted 
when the w-quark is boosted (giving rise to the finite- ir Bj 
component), while another component looks the same in 
all frames (contributing only at XBj =0). In 3 + 1- 
dimensional examples (similar to the one above) this 
means that the hadron does not become a "pancake" 
as p — > oo. Only part of the wave function is Lorentz 
contracted while some rest retains its spherical shape. 

There is nothing wrong or inconsistent about this re- 
sult. It is only another example where naive intuition 
(in this case based on the Lorentz transformation prop- 
erties in free field theory) fails. Interacting fields may 
have very complicated transformation properties which, 
in the above example, cause that part of the wavefunc- 
tion does not Lorentz contract in the same way as free 
fields do. This also "explains" why higher-twist distribu- 
tions are more vulnerable to these effects, since, in the 
parton model, they contain interaction terms explicitly. 

From the point of view of null plane quantization, the 
(5-functions are a manifestation of zero-modes 



4 30 



the hadronic wave function (see Appendix A). Higher- 
twist parton distributions are defined through correla- 
tions which involve bad currents (i.e. t/>( — ) compo- 
nent of the spinors). When solving the constraint equa- 
tion for tptj one has to specify boundary conditions at 
x~ = ±oo. Non-zero values at the boundary are of course 
related to (5-functions in the Fourier transform. Another 
(not necessarily independent) connection to the failure of 
canonical light-cone quantization (without zero modes) is 
observed by considering the singular part of e(x). In per- 
turbation theory one finds that the derivative of the non- 
covariant 7 + /p + term in the fermion self-energy with re- 
spect to the quark mass is proportional to the 5-function 
coefficient in e(x). Again this should not be a surprise, 
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since those terms are caused by an improper treatment of 
the zero modes in canonical light-cone perturbation 
theory [|l| — unless one uses several Pauli-Villars reg- 
ulators (more than in a manifestly covariant approach) 

EE- 

In practical applications, most sum rules for higher- 
twist distributions a re p robably useless since one expects 
them to diverge [ p3|j34[| . The discussion of e(x) in this 
work should therefore be considered only a pedagogical 
example. One exception seems to be gi- From a Regge 
pole analysis one expects the Burkhardt-Cottingham 
sum rule to converge H]. Also, as shown in Section EI], 



perturbation theory does not indicate any divergences in 
QCD. Furthermore, Regge pole analysis [Q as well as per- 
turbative considerations (Section IV) suggest validity of 
the sum rule. However, no strict predictions for gi have 
been made on a non-perturbative basis so far and at this 
point one has to wait for the experimental result to see 
whether the (72-sum rule is valid. 

For strange or heavier quarks, the sum rule seems 
to be most certainly violated. The perturbative analysis 
in this work only yielded a violation proportional to the 
quark mass and thus one might be tempted to assume 
validity of the Hl sum rule for u and d quarks. How- 
ever, it is not clear whether other effects, which have not 
been considered in this work, lead to a violation for light 
quarks as well. 

Even if these sum-rules are violated, the OPE remains 
valid. What needs to be added is a subtraction in the 
dispersion relation. Note that such a mechanism cannot 
be used to "explain" the spin crisis, since g\ is leading 
twist and thus a subtraction in the dispersion relation 
would be in conflict with unitarity constraints. 



APPENDIX A: E(X) ON A FINITE INTERVAL 

The 5-function at the origin in e(x) strongly sug- 
gests a relation to zero modes in light-cone quantization 
[§4],|Z|[30||. Due to the severe infrared singularities in 
light-cone quantization (the kinetic energy diverges as 
k + — > 0) it turns out to be necessary to go on a finite 
interval. Furthermore, following Lenz et al. [0], we will 
use a tilted coordinate system 



x = x 



L 



(Al) 



where the " refers to the "usual" light-cone coordinates 

(A2) 



(x^x 1 ) N2 



The notation in this Appendix will thus be different from 
the rest of the paper in order to allow for a direct com- 
parison with Ref. p4[ , which is strongly recommended to 
supplement this Appendix. 

This implies for example for scalar products 



a ■ b — a + 6_ + a_6 4 



(a_ + € -a + \ 
(A3) 



L refers to the length of the interval and we will im- 
pose periodic boundary conditions in x~ . The canoni- 
cally conjugate momentum p_ thus becomes quantized 



P-r 



2n 7r 



(A4) 



For the linear cr-model ( [5.l| ) one thus finds to order g 2 
ft 1 



e(k-) cx / dk+ 



(k 2 — m 
1 

dx 





2 + ief (p - k) 2 -m 2 + ie 
1 3 m 2 (l - x)" 



2£>3/2 4 £,5/2 



-1 -3/2' 



(A5) 



where 



D = — (k— — a;p_) 2 — x(l — x)m 2 + m 2 . (A6) 

Of course k- as well as p- take on only discrete val- 
ues (A. 4). Here c is some numerical constant propor- 
tional to g 2 . One immediately recognizes the similarity 
to Eq. (3.8), which is not completely accidental, since the 
infinite momentum boost also gives some kind of regu- 
larized representation for e{x). The canonical light-cone 
limit for e(fc_) is obtained by taking eLm 2 — > while 
keeping and n p fixed (fc_ = 2-KTik/L, p_ = 2irn p /L). 
Using 



dx f(x) 



(A7) 



this yields (y 



k- 



e(fc_) 



V2- 



1 1 
y z m z 



[y(l — y)m 2 — m 2 ] 
for rik 7^ 0. However, for iik = 0, 

e(0) - ' ' 



T 8(y)0(l-y) , (A8) 



(A9) 



i.e. the zero mode contribution diverges. 

The mere fact that e(0) diverges reflects the existence 
of the 5-function in the continuum calculation. However, 
if one actually wants to make a quantitative comparison 
with Section III, one has to perform a more sophisticated 
continuum limit where also the hadron momentum n^. (in 
discrete units) approaches infinity. 
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